Matter-wave bright solitons of "^Li gas in an expulsive potential 
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A Lagrangian based method is used to derive an analytical model for studying the dynamics of 
matter-wave bright soliton created in a harmonic potential which is attractive in the transverse 
direction and expulsive in the longitudinal direction. By means of sech trial functions and a Ritz 
optimization procedure, evolution eqautions are constructed for width, amplitude and nonlinear 
frequency chirp of the propagating soliton of the atomic condensate. Our eqaution for the width is 
an exact agreement with that of Carr and Castin [Pftj/s.iiew.j4 66, 063602(2002)] , obtained by more 
detailed analysis. In agreement with the experiment of Paris group [Science 296, 1290(2002)], the 
expulsive potential is found to cause the soliton to explode for A^|as| = 0.98, A'^ being the number 
of atoms in the condensate and a^, the scattering length of the atom-atom interaction. 

PACS numbers: 03.75.Fi, 05.45.Yv 



Solitons are localized waves that can propagate over 
long distances with neither attenuation nor change in 
shape. Such waves are formed only when their dispersion 
is compensated by the nonhnear effects of the medium. 
Sohtons appear in many diverse physical situations in- 
cluding waves in shallow water , pulse propagation in 
optical fibers and plasma waves. Since the experimen- 
tal observation of Bosc-Einstein condensates (BEG) in 
1995, it was realized that we can also have matter-wave 
solitons P]. In this case, the nonlinearity is produced 
by binary atom-atom interactions leading to the mean 



field U{r) 



where is the s-wave scattering 



length for the atom-atom scattering and m, the atomic 
mass. Here ^j{r,t) is the wavefunction or order param- 
eter of the condensate. For > we get dark solitons 
while for < we have bright solitons. A dark soliton 
is a notch on the BEG with a characteristic phase steep 
across it. On the other hand, bright soliton is a peak. 

Almost simultaneously two different groups, one at 
Ecole Normale Superiore in Paris and other at Rice 
University produced bright solitons from a trapped 
^Li BEG in the internal atomic state \F = l,mF = 1) 
by continuously tuning the scattering length from a pos- 
itive to a negative value with the help of Feshbach res- 
onance iJl] induced by an applied magnetic field. In the 
Paris experiment, the spherical trapping of the BEG was 
adiabatically deformed into a cylindrical geometry and 
the condensate was finally released into a horizontal ID 
waveguide such that the resulting force on the atoms can 
be conveniently represented by —muj'^z, where ujz is the 
imaginary frequency along the z direction. Understand- 
ably, this imaginary frequency is due to an offset mag- 
netic field that produces an expulsive harmonic potential 
— ^TOW^z^. The main effect of this expulsive term in the 
evolution equation is that the center of mass of the BEG 
accelerates along the longitudinal direction. Thus within 
the framework of mean-field approximation the Paris ex- 



periment can be modelled by a 3D Gross-Pitaevskii equa- 
tion 



t/- = 



with 



(1) 



(2) 



the confining potential. Here u;± is the frequency of the 
radial trapping and the wave function ili{r,t) is normal- 
ized to the number of particles N in the condensate such 
that 



I if; |2 dr = N. 



(3) 



In the case strong cylindric radial confinement it will be 
convenient to work with a quasi one-dimensional (QID) 
form of (1) and thereby study the dynamics of matter- 
wave bright solitons created in a harmonic potential 
which is attractive in the transverse direction but ex- 
pulsive in the longitudinal direction. We shall envisage 
such a study with particular emphasis on the role of ex- 
pulsive potential in causing explosion as observed in the 
Paris experiment 2]- 

There exists a number of detailed studies in respect 
of this. Remarkable among them are the works of Garr 
and Gastin Jj and of Salasnich 6]. Our object in this 
work is to demonstrate how a simple analytical model 
could be used to physically realize the problem of soli- 
ton explosion. We shall show that the critical number 
of atoms in the soliton just before explosion is in exact 
agreement with the numbers predicted by the sophisti- 
cated works in refs. 5 and 6 . To that end we first obtain 
QID equation from (1). We write this equation in terms 
of dimensionless variables defined by 



ij{r,z,t) = u{p,s,T)/a§ 
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with ao = 



-, the size of the ground state solution 



2 



of the noninteracting GP equation. This gives 
iur + \v^u -lip-"- \y) u-^\u\^u = 0, A, 



2 2 
It is obvious that 



= N. 



Using a separable ansatz 

u{p,s,t) = (t){p)^{s,T) 
(4a) can be written in the form 



(4a) 



(4b) 



(5) 



1 



1 



? V 2 2 J aa 



1 



(6) 



where stands for the Laplacian in the radial coordi- 
nate. In (6) the subscripts on ^ stands for partial deriva- 
tive with respect to that particular independent variable. 

More specifically, ^2s = q^- This equation shows that 
the presence of atom-atom interaction does not permit 

clcarcut separation variables. However, the fourth term 
in equation (6) is quite small such that (j) may be assumed 
to satisfy 



1-72^ , 1 2, 



(7) 



with LUp related to lu±_ by a scale factor. Equation (7) 
represents the well-known eigenvalue problem for the two 
dimensional harmonic oscillator with the ground sate so- 
lution given by 

Thus (6) can be written in the form 

+ + ^xy^ - ^\em = (8) 

z z ao 

The low-frequency vibration along the z direction is quite 
unlikely to excite the two dimensional bosonic oscillator 
from its ground state. In view of this we multiply (8) by 
(f)(p* and integrate over the p coordinate to get 



2„2. 2.a.|^|,^^^^^^ 



ao 



Equation (9) can be written in a more convenient form 
by using the change of variable 



^(s,T)=x(s,T)e-*"''". 
Using (10) in (9) we get 



(10) 



iXr + lx2s + Ixyx - ^\X\\ = (11a) 
z z ao 



with 



/ + 00 
I 
— oo 



\xfds = N/tt. 



(lib) 



The initial-boundary value problem in (11a) can be 
converted to a variational problem in a rather straight- 
forward manner. For example, it is easy to see that the 
action principle 

S J J C{x,X*,Xs,x:,Xr,X*)dsdT (12) 

is equivalent to (11a) with the Lagrangian density 

^ = ^ ixxr - x*Xt) - ^xyxx* + —x\*^ + \x*sXs- 

z z ao z 

(13) 

As a solution for (11a) we introduce the trial function 
X{s, t) = A(T)scch (,s/a(r)) e^^(^)^' . (14) 

with a complex amplitude A{t). Here a(T) stands for the 
width of the distribution and 6(t), the frequency chirp. 
In terms of trial function in (14) we obtain a reduced 
variational problem. 



i j (£) dT = 







with 



/'+CC 
-oo 



Cs ds. 



(15) 



(16) 



Here Cs stands for the values of C when (14) is substi- 
tuted in (13). We perform the integration in (16) to get 



. , 7r2 



C^ia {AA* - A* At) + —AA*bta^ Xja^AA* 



6 



12 



—^aA^A*^ + — — + —b^a^AA*. (17) 
3ao 3a 3 ^ ' 

The vanishing conditions of ^3a^ ^^"^ ^W' 

yield 



2 

6 12' 

,2„3 A* 



2iaA* + iatA* + ^a^btA* - ^Xla^A* 



87ra,, . ,^^2 A* tt"^ 



aAA*' + — + —b^a'^A* = , (18a) 
3ao 3a 3 ' v / 



2iaAt + iatA - —a^btA+ —Xla^A - 
6 12 

^aA*A' - - - -b'a'A = , (18b) 
3ao 3a 3 ^ \ J 



i {AA* - A* At) + —bta^AA* - —XyAA* ^ 



12 



3ao 



'A^A*^- 



AA* 
"3^ 



TT^Va^AA* = , (18c) 



3 



and 



dt 



a^AA") -Aba^'AA* = 



(18d) 



Equations (18a) and (186) can be combined to get 

a|Ap = Constant . (19) 

From the normalization of the trial function, the constant 
in (19) can be identified with N/2tt such that 



a\A[ 



N 
2^ 



(20) 



Equations (18a) - (18c) and (20) can be combined to 
write 

a%r-T:XW-^ ^ + 20^52 ^0. 21 

2 TT^ ao a TT^a-^ 



From (18d) and (20) we have 

b 1 '^'^ 
~ 2a eft ■ 



(22) 



Using (22) in (21) we obtain a second-order differential 
equation for a(T) given by 



a2T — A^a 



4 Na, 1 



TT^ gq a?- TT^a^ 



0. 



Equation (23) can be integrated to get 
1 



with 



-a^ + Y(a) = E 



2 a TT^ a^ 



(23) 



(24) 



(25) 



Here is a constant of integration and P 



4 Na^ 



The extrema of the potential function in (25) are de- 

V{a 
da 

TT^Xla'^ - 7r2|P|a + 4 = . (26) 



termined from '^'^Jf' = 0. This gives 



It is of interest to note that (26) is an exact agreement 
with equation (19) of Carr and Castin who studied the 
dynamics of matter-wave bright soliton with special at- 
tention on the variation of Gross-Pitaevskii energy func- 
tional. However, as opposed to the result in ref. 5, the 
width a in (26) is a function of time. This implies that 
(26) holds good for any value of t. Whether these ex- 
trema are saddle points, minima, maxima is determined 
by 



d^Vja) 
da? 



2\P\ 



12 



(27) 



In writing (26) and (27) we have introduced P — —\P\ to 
carry out the subsequent analysis by using only numerical 
values of a^. We note that for bright solitons the nonlin- 
ear term is attractive and the scattering length a^ < 0. 
The critical point for collapse caused by the mean field, 
in the case < 0, for explosion due to expulsive poten- 
tial may be obtained by simultaneously solving (26) and 

= 0. We have 



(27) for a and \P\ when 



37r2A? 



and \P\ = 



16 



(28) 



Equation in (28) can be combined to get 

iV|a,| = (4/3)'/^7ri/2Ai/2ao . (29) 

Using the experimental parameters ao = 1.43/xm, Az = 
^ we have N\ai\ — 0.98 for which the expulsive poten- 
tial causes the soliton to explode axially. This theoretical 
result is in agreement with the prediction of Khaykovich 
et al to within 10%. We, therefore, conclude by noting 
that there are distinct advantages to viewing nonlinear 
dynamics of bright solitons within the framework of a 
Lagrangian based approach. 
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